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Evaluation  of  Habitat  Suitability  Index  Models  for  Riverine 
Life  Stages  of  American  Shad,  with  Proposed  Models 
for  Premigratory  Juveniles 


by 


Robert  M.  Ross,  Thomas  W.  H.  Backman,^  and  Randy  M.  Bennett 


U,S.  Fish  and  Wildlife  Service 
National  Fisheries  Research  Center — Leetown 
Natwnal  Fishery  Research  and  Development  Laboratory 
Rural  Delivery  #4,  Box  63 
Wellsboro,  Pennsylvania  1 6901 


Abstract.  Field  evaluations  of  existing  habitat  suitability  index  (HSI)  models  for 
spawning  adults,  eggs,  and  larvae  of  American  shad  (Alosa  sapidissima)  were  conducted 
in  1990-92;  initial  models  for  juveniles  in  nursery  habitats  were  developed.  Fish 
abundance  in  various  habitats  of  the  upper  Delaware  River  was  quantified  by  (1) 
observation  of  adult  spawning  activity,  (2)  collection  of  eggs  and  larvae  with  metered 
plankton  and  drift  nets,  and  (3)  enumeration  of  juveniles  by  underwater  observation  and 
seining  techniques.  Regression  analysis,  principal  component  analysis,  and  range 
analysis  were  used  to  relate  abundance  to  an  array  of  physical  habitat  variables 
potentially  influencing  fish  distributions. 

The  mean  surface  water  temperature  model  for  spawning  adults  (maximum  suitability 
of  14-20®  C)  was  supported  at  its  lower  limit,  but  the  data  suggest  an  upper-limit 
extension  to  24.5®  C.  The  mean  water  velocity  model  (maximum  suitability  of  0.3-0.9  m/  s) 
was  not  supported,  particularly  at  its  lower  limit  (0-0.3  m/s),  where  adults  spawned 
with  equal  frequency.  Maximum  suitability  at  the  upper  limit  did  not  extend  beyond 
0.7  m/  s.  Spawning  activity  was  observed  in  all  habitat  types  examined  but  was  greatest 
in  runs  and  lowest  in  pools  and  riffle  pools. 

The  HSI  model  for  eggs  and  larvae  in  relation  to  surface  water  temperature  (maximum 
suitability  of  15-25®  C)  was  supported  for  eggs  and  prolarvae.  Surface  water  temperature 
was  not  predictive  of  postlarval  relative  density,  however,  because  no  reduction  in  density 
occurred  at  the  upper  thermal  limit  (26-27®  The  temporal  variables  date  and  time  of 
day  influenced  egg  and  larval  densities.  Prolarvae  were  also  afrected  positively  by  sample 
depth  (among  five  remaining  physical  habitat  variables),  while  postlarval  densities  were 
influenced  positively  by  water  depth  and  temperature.  Only  postlarvae  showed 
differential  habitat  use:  riffle  pools  were  selected  over  riffles,  channels,  pools,  and  slopes. 

No  HSI  model  was  previously  developed  for  juvenile  American  shad  in  riverine  habitats. 
We  found  four  physical  habitat  variables  that  were  correlated  with  juvenile  abundance: 
water  temperature,  dissolved  oxygen  (covariates),  river  depth,  and  turbidity.  The 
influence  of  these  variables  was  habitat-specific  (six  types).  Preliminary  HSI  models  were 
developed  for  three  of  these  variables.  There  was  no  overall  effect  of  habitat  type  on 
juvenile  relative  abundance,  indicating  that  juveniles  use  a  wide  variety  of  habitat  types 
to  their  advantage  in  many  nursery  areas.  In  contrast  to  earlier  life  stages  and  spawning 
adults,  premigratory  juveniles  appear  to  be  habitat  generalists. 


^  Present  address:  Columbia  River  Inter-Tribal  Fish  Commission, 
729  Oregon  Street,  Suite  200,  Portland,  Oreg.  97232. 
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Suitability  ranges  for  some  physical  variables  were  narrower  or  broader  than  previously 
assumed;  moreover,  the  influence  of  these  variables  was  found  to  be  habitat  specific. 
Therefore,  American  shad  HSI  models  may  be  applied,  with  caution,  on  a  comparative 
basis  among  river  systems.  Knowledge  of  the  range  and  proportion  of  habitat  types  in 
each  river  system  considered  would  be  helpful  in  applying  the  models. 

Key  words:  Alosa  sapidissima,  American  shad,  habitat  suitability  models,  larval  fishes, 
reproduction,  river  ecology. 


Habitat  suitability  index  (HSI)  models  were  de¬ 
veloped  in  the  1970's  and  1980's  to  facilitate  impact 
assessment,  project  planning,  habitat  manage¬ 
ment,  and  the  understanding  of  species-habitat 
relations  (U.S.  Fish  and  Wildlife  Service  1980; 
Schamberger  et  al.  1982;  Stier  and  Crance  1985). 
These  models  provide  information  on  species-habi¬ 
tat  correlations  that  may  be  strengthened  with 
further  study.  Many  HSI  models  were  developed  on 
the  basis  of  known  literature  and  the  Delphi  proc¬ 
ess  (Crance  1987).  The  Delphi  process  relies  on 
extensive  experience  with  the  concerned  species  by 
the  individuals  queried.  Habitat  field  studies  were 
not  directly  employed  to  generate  a  data  base  from 
which  many  HSI  models  were  developed.  Such 
models  are  known  as  category-one  suitability  index 
curves  (Bovee  1986). 

The  HSI  models  for  American  shad  (Alosa 
sapidissima)  were  developed  for  spawning  adults, 
eggs,  and  larvae  in  riverine  habitats,  but  models  for 
juveniles  were  limited  to  estuarine  life  during  win¬ 
ter  and  spring  (Stier  and  Crance  1985).  Habitat 
variables  deemed  important  to  American  shad  dur¬ 
ing  riverine  components  of  their  life  history  were 
mean  surface  water  temperature  during  spawning 
season  or  egg-larval  development  period  (spawn¬ 
ing  adults,  eggs,  and  larvae)  and  mean  water  veloc¬ 
ity  during  spawning  season  (spavraing  adults). 
Other  variables  such  as  depth,  substrate,  and  cover 
were  thought  not  to  appreciably  influence  habitat 
use  by  American  shad. 

The  purpose  of  our  study  was  to  evaluate  the 
existing  HSI  models  for  American  shad  in  riverine 
habitats  by  generating  suitability  index  (SI)  curves 
for  the  Delaware  River,  one  of  the  most  productive 
American  shad  river  systems  of  the  eastern  United 
States.  We  compared  our  curves  with  existing  SI 
curves  (Stier  and  Crance  1985)  and  evaluated  the 
models  with  statistical  tests  of  correlation  to  habi¬ 
tat  variables.  We  also  developed  habitat  utilization 
relations  for  juvenile  (premigratory)  American 
shad  in  nursery  habitats,  and  we  propose  HSI  mod¬ 
els  for  this  previously  ignored  life  stage.  The  cate¬ 


gory-two  suitability  index  curves  (Bovee  1986)  gen¬ 
erated  by  this  study  are  considered  to  have  broader 
applicability  than  those  of  category  one,  but  they 
are  not  habitat  preference  curves. 

Methods 

Study  Area 

Field  work  was  conducted  on  the  Delaware 
River  between  the  Delaware  Water  Gap  (river 
kilometer  341,  neeir  East  Stroudsburg,  Pennsylva¬ 
nia)  and  Callicoon,  New  York  (river  kilometer  488; 
Fig.  1).  The  lower  reaches  of  this  river  segment 
fall  within  the  boundaries  of  the  Delaware  Water 
Gap  National  Recreational  Area,  while  upper 
reaches  fall  within  the  Upper  Delaware  Scenic 
and  Recreational  River,  both  administered  by  the 
National  Park  Service  (NPS).  Fisheries  jurisdic¬ 
tion  includes  the  NPS  and  states  of  Pennsylvania, 
New  Jersey,  and  New  York. 

Spawning  Adults 

We  quantified  spawning  of  adult  American 
shad  by  counting  spawning  splashes  over  5-min 
intervals  between  darkfall  and  0100  h  (Eastern 
Daylight  Saving  Time)  from  16  April  to  20  June 
1991  and  13  April  to  9  July  1992.  To  verify  that 
such  splashes  were  made  by  spawning  American 
shad,  we  used  spotlights  from  an  anchored  boat  at 
two  routinely  observed  spawning  sites.  Over  a  1-h 
period  at  one  site  and  0.5  h  at  another,  spotlights 
were  activated  for  each  splash  heard  within  about 
5  m  of  the  boat.  The  number  and  species  of  fish 
observed  and  their  behavior  were  recorded.  We 
assume  that,  though  every  splash  may  not  repre¬ 
sent  actual  spawning  and  every  spawning  may  not 
be  accompanied  by  a  splash,  the  level  of  surface 
activity  is  strongly  correlated  with  actual  spawn¬ 
ing.  Eleven  sites  (Table  1)  representing  five  habi¬ 
tat  types  (Table  2)  were  sampled  1  to  17  times 
each  by  stationing  two  listeners  on  shore  or  in  an 
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anchored  boat  at  each  site.  The  average  of  the  two  5  IB),  eind  turbidity.  Turbidity  was  obtained  from 
counts  of  all  splashes  heard  from  a  single  location  a  turbidity  meter  reading  (Markson  nephelome- 
was  recorded.  After  each  count  the  following  habi-  ter)  of  a  refrigerated  25-mL  surface  water  sample 
tat  variables  were  recorded:  river  depth  (10  m  within  24  h.  Habitat  type  was  also  determined 
from  shore  or  at  boat  anchorage),  middepth  cur-  and  recorded  at  each  count  location, 
rent  velocity  10  m  from  shore  or  at  anchorage  We  used  principal  component  aneilysis  (PCA; 
(Marsh  McBimey  Model  201D),  surface  water  SAS  1987)  to  explore  relations  between  spawning, 
temperature,  dissolved  oxygen  (DO;  YSI  Model  five  physical  habitat  variables,  and  five  habitat 


Table  1.  Access  points  and  life  stages  sampled  for  American  shad  (Alosa  sapidissima)  on  the  upper 

Delaware  River,  1990-92. 


Site  access 

River  kilometer 

Spawning  adults 

Life  sta^re 
Eggs 

sEimpled 

Larvae 

Juveniles 

Kittatinny  (New  Jersey) 

341 

X 

Worthington  State  Forest  (New  Jersey) 

346 

X 

Smithfield  Beach  (Pennsylvania) 

351 

X 

X 

X 

X 

Pc>xono  (New  Jersey) 

354 

X 

Dingman's  Ferry  (Pennsylvania) 

383 

X 

X 

X 

X 

Milford  Beach  (Pennsylvania) 

397 

X 

X 

X 

X 

Mongaup  (New  York) 

418 

X 

Zane  Grey  (Pennsylvania) 

447 

X 

X 

X 

X 

Narrowsburg  (Pennsylvania) 

467 

X 

Damascus  (Pennsylvania) 

479 

X 

Callicoon  (Pennsylvania) 

488 

X 
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types.  Observations  were  grouped  into  three 
classes  by  number  of  splashes:  0,  1-25,  and  >25. 
The  average  value  for  principal  components  (PC) 
1  and  2  of  each  group  was  plotted  on  the  PC1-PC2 
grid,  along  with  vectors  for  each  habitat  variable. 
Some  degree  of  correlation  among  independent 
variables  is  assumed  in  PCA,  but  not  in  regression 
analysis  (Lindeman  et  al.  1980).  To  further  exam¬ 
ine  relations  with  specific  habitat  variables,  we 
used  linear  and  quadratic  regressions  (SAS  1987) 
with  mean  number  of  spawning  splashes  as  the 
dependent  variable  and  the  above  physical  factors 
as  independent  variables.  A  null-hypothesis  rejec¬ 
tion  level  of  P  <  0.05  was  set  for  these  and  all 
subsequent  statistical  tests.  We  used  the  Duncan 


multiple  range  test  (SAS  1987)  to  rank  and  deter¬ 
mine  significant  differences  in  spawning  activity 
among  five  habitat  types. 

Eggs  and  Larvae 

We  sampled  American  shad  eggs  and  larvae  at 
four  locations  (Table  1),  representing  six  habitat 
types  (Table  2),  with  calibrated  plankton  nets  in 
1990  (8  May  to  8  July),  1991  (16  April  to  20  Jxme), 
and  1992  (13  April  to  9  July),  from  1000  to  0230  h. 
Four  types  of  plankton  nets  (all  500-pm  mesh  size) 
were  used:  (1)  61  x  25-cm  (diameter)  drift  nets,  (2) 
35  X  35  X  61-cm  benthic  sled  (inverted  "U"  opening) 
with  1-L  bucket,  (3)  2.5  x  0.5-m  (diameter)  bongo 


Table  2.  Description  and  distribution  of  nine  habitat  types  among  the  four  primary  oecess  sites  for 
American  shaders,  larvae,  and  juveniles  on  the  upper  Delaware  River.  Five  of  these  habitats  (channel, 
shallows  with  submerged  aquatic  vegetation  [SAl^,  run,  riffle  pool,  and  pool)  were  also  assessed  for 
spawning  adults. 

Sites  where  sampled^ 

Habitat  type  Eggs  and  larvae  Juveniles  Description  of  habitat 


Channel 

S,D,M 

D 

Eddy 

D 

Slope 

S 

SAV  shallows 

S,D,M 

S,D,M, 

Run 

Riffle 

S,M,Z 

S,M,Z 

Riffle  pool 

Z 

S,M,Z 

Pool 

Z 

Tributary 

W,M 

Deeper  (2-5  m)  portion  of  river  where  greatest  current 
velocities  found;  typically  the  middle  third  of  river 
width;  little  if  any  macrophytic  vegetation  anchored 
in  substrate 

Zone  of  reversed  current  direction  lateral  to  channel; 
includes  shore  areas  with  SAV  but  narrow  shallow 
zone  (rapid  drop  to  >2  m) 

Segment  of  river  width  shoreward  of  channel,  1-3.5  m 
deep,  but  not  including  shallow  inshore  zone;  some 
attached  macrophytes  may  be  present,  but  not  at  high 
densities  typically  found  inshore 

Inshore  zones,  0-1.5  m  deep,  where  high  densities  of 
anchored  macrophytic  vegetation  and  low  current 
velocities  typically  occur;  sampled  quantitatively 
only  for  juveniles 

Midriver  stretch  (spawning  adults  only)  of  relatively 
shallow  (0.5- 1.5  m)  water  of  moderate  to  high  current 
velocity  (0.3-0. 7  m/  s) 

Zone  of  shallow  (<0.5  m)  swift  (typically  >0.5  m/  s) 
water  with  surface  turbulence  and  aeration;  SAy  if 
present,  typically  short  (<0.25  m) 

Zone  immediately  downstream  of  a  riffle,  where  water 
deepens  suddenly  and  currents  become  variable  in 
velocity  and  direction 

Relatively  large  river  segment,  where  river  width 
bulges  and  current  velocities  drop  (<0.25  m/  s);  typically 
just  upstream  of  a  series  of  riffles 

Mouth  and  delta  of  a  tributary  stream,  characterized  by 
relatively  cold  water  temperatures,  steep  water 
temperature  gradients,  and  visibly  dense  suspended 
material  ("drift") 


S  =  Smithfield  Beach,  D  =  Dingman's  Ferry,  M  =  Milford  Beach,  W  =  Delaware  Water  Gap,  and  Z  -  Zane  Grey. 
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net  with  1-L  buckets,  and  (4)  6.1  x  2.4-m  plankton 
seine  with  2.4  x  1.2  x  2.4‘m  bag  and  1-L  bucket. 
Immediately  before  or  after  (or  during  in  the  case 
of  drift-net  and  plankton-seine  sets),  the  following 
physical  variables  were  measvired  and  recorded: 
average  river  depth,  sample  depth  (from  river 
depth  and  gear  used),  water  temperature,  DO,  cur¬ 
rent  velocity,  and  turbidity 

Drift  nets  were  deployed  in  triplicate  in  riffle 
zones  only  and  were  secured  with  horizontal  rein¬ 
forcing  bars  attached  to  heavy  metal  stakes.  Veloci¬ 
ties  through  nets  (used  to  calculate  volumes  of 
water  sampled)  were  measured  with  a  Marsh 
McBimey  Model  201D  current  meter  at  the  begin¬ 
ning  and  end  of  each  1.1-3.7-h  set  and  averaged. 
The  benthic  sled  was  equipped  with  a  propeller- 
style  digital  flow  meter  calibrated  with  the  above- 
described  current  meter.  It  was  deployed  from  a 
boat  with  a  1:1  winch  mounted  to  a  chest-high  davit 
and  15  m  of  3-mm  steel  rope.  The  sled  was  towed 
at  a  speed  that  allowed  it  to  skid  lightly  on  the  river 
bottom.  The  bongo  plankton  seine  was  outfitted 
similarly  but  towed  at  a  constant  speed  just  below 
the  surface. 

The  large  plankton  seine  was  also  outfitted  with 
a  calibrated  flow  meter  and  held  open  by  attach¬ 
ment  to  a  horizontal  3.7-m  steel  rod.  The  preferred 
method  of  deployment  was  about  8  m  behind  the 
stem  of  the  boat,  which  had  to  be  anchored  firmly 
to  a  large  rock  or  similar  object.  Firm  anchorage 
could  not  be  fovind  at  one  site,  however,  so  dual 
grapnels  were  deployed  off  the  bow  to  reduce  boat 
drift  and  to  increase  flow  through  the  net  while 
dragging. 

Volumes  of  water  sampled  per  set  by  these 
techniques  ranged  as  follows:  drift  net  20-964  m^, 
benthic  sled  1-93  m^  bongo  net  51-316  m^,  and 
plankton  seine  11-4,806  m^.  All  samples  were 
carefully  removed  from  nets  or  buckets  and  fixed 
in  10%  formalin  for  transport  to  the  laboratory, 
where  they  were  transferred  to  70%  ethanol  for 
preservation.  All  nonadhesive  eggs  2.3-3.8  mm 
total  diameter  (after  immersion  >10  min  in  deion¬ 
ized  water)  with  nuclear  diameter  <66%  of  total 
diameter  were  identified  as  American  shad  eggs. 
These  criteria  were  necessary  to  distinguish  Ameri¬ 
can  shad  eggs  from  those  of  the  white  sucker 
(Catostomus  commersoni),  which  spawns  at  the 
same  time  but  produces  fertile  eggs  with  a  differ¬ 
ent  frequency  distribution  of  egg  to  nucleus  size 
(R.  M.  Ross  and  Bennett  1993).  American  shad 
yolk-sac  (prolarvae)  and  post-yolk-sac  (postlar¬ 
vae)  larvae  were  distinguished  from  anatomically 


similar  larvae  and  identified  by  number  of  preanal 
myomeres,  number  of  postanal  myomeres,  rela¬ 
tive  preanal  length,  and  pigmentation  patterns 
(Lippson  and  Moran  1974;  Fuiman  1979;  Holland- 
Bartels  et  al.  1990).  We  captured  only  larvae 
<20  mm  in  total  length. 

Egg  and  larval  densities  were  calculated  by 
dividing  the  number  of  eggs  or  larvae  caught  by 
the  volume  of  water  sampled  by  the  net.  Because 
four  different  sampling  techniques  were  used  to 
obtain  all  three  life  stages,  the  possibility  of  dif¬ 
ferent  catchability  of  the  more  mobile  postlarval 
stage  was  considered.  We  assume  no  such  differ¬ 
ence  for  eggs  and  the  weak-swimming  prolarvae. 
We  calculated  sampling  effort  by  gear  type  and 
habitat  type  (Table  3).  Differential  catchability  of 
postlarvae  was  evaluated  indirectly  through  com¬ 
parison  of  observed  postlarval  densities  in  the 
same  habitat  for  different  gear  types.  Two  com¬ 
parisons  met  these  requirements:  (1)  plankton 
seine  and  bongo  net  in  near-surface  channel  habi¬ 
tat  and  (2)  drift  net  with  bongo  net  in  riffle  habitat. 
Postlarval  densities  in  various  habitats  were  also 
compared  by  gear  type  to  determine  whether  dif¬ 
ferent  conclusions  might  be  reached  depending  on 
sampling  gear  employed. 

Relations  between  egg  and  larval  density,  physi- 
cal  habitat  variables,  and  habitat  type  were  ex¬ 
plored  with  PCA  by  superimposing  habitat  type 
and  life  stage  on  the  PC  1-2  xmiverse  generated  by 
physical  habitat  variables  (available  habitat).  Life- 
stage  density  values  were  weighted  by  presence  or 
absence  (1  or  0)  criteria  to  achieve  better  separation 
of  life  stages.  Density  values  were  also  used  as 
dependent  variables  in  linear  and  quadratic  regres¬ 
sion  analyses  (SAS  1987)  with  the  following  con¬ 
tinuous  habitat  variables:  sample  depth,  river 
depth  where  sample  was  taken,  water  tempera¬ 
ture,  DO,  current  velocity,  and  turbidity.  Also,  the 
effect  of  habitat  type  on  egg  and  larval  abundance 
was  tested  by  subjecting  egg  and  larval  densities  to 
analysis  of  variance  (ANOVA;  SAS  1987)  with  four 
habitat  types  as  class  variables  (Table  2;  heavily 
vegetated  shallows  excluded  because  metered  nets 
were  not  employed  there,  riffle  pool  excluded  be¬ 
cause  of  small  sample  size).  We  used  the  Duncan 
range  test  (SAS  1987)  to  rank  mean  den¬ 
sities  for  each  habitat  type.  Scatter  plots  were  gen¬ 
erated  to  illustrate  egg  and  larval  densities 
(nonzero  values  only)  as  a  function  of  each  meas¬ 
ured  physical  habitat  or  temporal  variable. 

Inshore  habitat  zones  with  high  densities  of  an¬ 
chored  macrophytic  vegetation  (submerged  aquatic 
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vegetation,  or  SAM  excluding  algal  mats)  were  not 
sampled  by  metered  plankton  nets.  Rather,  this 
habitat  type  was  sampled  in  1990  with  the  6.1-  x 
2.4-m  plankton  seine  (unmetered)  by  walking  the 
net  from  shore  to  a  maximum  depth  of  1.5  m  and 
hauling  it  in  to  shore  for  egg  and  larval  collection. 
Water  volumes  sampled  varied  considerably  with 
each  haul,  but  20  hauls  were  made  in  shallow  SAV 
habitats  at  three  different  sites  between  22  May 
and  8  July  1990. 

Juveniles 

Young-of-year  juvenile  American  shad  were 
quantified  in  six  habitat  types  at  five  access  sites 
with  two  different  methods:  seine  hauls  and  quan¬ 
titative  underwater  observations  (Table  2).  For 
seine-haul  determinations,  a  30.5-  x  2.4-m  (8-mm 
mesh)  blackened  beach  seine,  which  captured  only 
juveniles  >35  mm  total  length,  was  deployed  from 
the  stem  of  a  boat  or,  in  shsJlow  habitats,  by  walk¬ 
ing  the  seine.  Initial  sets  approximated  a  10-m 
square  before  hauling  toward  a  single  landing 
site  (on  shore  or  in  shallow  riffles  if  required); 
however,  in  habitats  subject  to  strong  currents 
(channels,  eddies,  and  riffle-eddy  zones),  square 
sets  quickly  became  linear  as  the  seine  was 
hauled  in. 

Gear  efficiency  was  estimated  in  SAV  shallow 
habitats  by  deploying  the  seine  by  foot,  as  described 
above.  A  duplicate  seine  was  then  deployed  just 
outside  the  first  as  a  block  net  and  the  inner  seine 
hauled  to  shore.  The  inner  net  was  reset  and  pulled 
in  two  more  times  before  the  outer  net  was  hauled 
to  shore.  Catchability  was  calculated  as  the  number 
of  juveniles  in  the  first  haul  divided  by  the  total 
number  offish  captured.  Catchability  studies  were 
conducted  four  times  in  SAV  shallow  habitats,  but 
they  were  not  feasible  in  other  habitat  types  be¬ 
cause  of  current-velocity  and  river-depth  limita¬ 
tions. 


Estimates  of  juvenile  relative  abimdance  were 
also  obtained  by  dive  counts  on  established  5-m- 
square  grids  in  each  habitat  type  at  all  access  sites 
(Table  3).  The  comers  of  all  grids  were  marked  by 
0.3-m  spikes  driven  into  the  substrate  and  flagged 
with  yellow  plastic  tape.  Before  any  seining  activity 
in  a  study  area,  counts  of  juveniles  were  made  by  a 
diver  with  mask,  fins,  snorkel,  and  slate.  Coimts 
were  made  by  swimming  up  to  the  perimeter  of  a 
grid  and  then  over  the  grid  as  needed.  Three  con¬ 
secutive  counts  were  used  to  determine  a  mean 
coimt  for  the  grid.  The  above  procedure  was  re¬ 
peated  after  10  min  to  assess  the  statistical  vari¬ 
ability  of  this  count  method. 

The  following  habitat  variables  were  obtedned 
from  the  sample  area  either  immediately  before  or 
after  (allowing  visibly  disturbed  water,  if  any,  to 
pass  first)  both  seine  hauls  and  dive-count  observa¬ 
tions:  average  river  depth,  water  temperature,  DO, 
average  current  velocity,  turbidity,  and  pervent 
SAV  cover.  Percent  SAV  cover  was  estimated  visu¬ 
ally  to  the  nearest  5%  in  each  dive-count  grid  or 
seine-haul  area. 

Seine-haul  and  dive-coimt  data  were  compared 
statistically,  both  overall  and  within  each  habitat 
type,  with  Pearson  correlation  coefficients  (SAS 
1987).  We  used  PCA  to  relate  juvenile  relative 
abundance  to  habitat  by  superimposing  habitat 
type  and  dive  or  seine  juvenile  coimts  on  the  uni¬ 
verse  of  vectors  representing  available  physical 
habitat.  The  relation  between  juvenile  abundance 
(dependent  variable)  and  each  measured  habitat 
variable  was  then  examined  by  regression  analysis 
(SAS  1987),  both  overall  and  by  habitat  type.  Both 
of  these  analyses  were  performed  independently 
with  dive-count  data  and  combined  dive-and-seine 
data  sets.  The  effect  of  habitat  type  on  juvenile 
relative  abundance  was  determined  by  ANOVA 
with  six  habitat  types  as  class  variables. 


Table  3.  Sampling  effort  in  the  upper  Delaware  River  for  American  shad  eggs  and  larvae  by  gear  type 
for  each  habitat  type.  Numbers  represent  volumes  of  river  water  sampled  in  cubic  meters,  with 
percentage  of  total  volume  for  each  habitat  type  (gear  type  under  "all  habitats")  in  parentheses. 


Gear  type 

Pool 

Channel 

Riffle^ 

Slope 

All  habitats 

Bongo  net 
Pleinkton  seine 

Drift  net 

Benthic  sled 

3,029  (90) 

347  (10) 

10,651  (45) 

12,115  (51) 

890  (4) 

3,972  (7) 

28,535  (49) 

25,709  (44) 

50  (0.1) 

290  (72) 

111  (28) 

17,942  (21) 

40,650  (47) 

25,709  (30) 

1,398  (2) 

Includes  riffle  pool  habitat. 
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Results 

Spawning  Adults 

Over  a  1-h  observation  period  at  one  site  and 
0,5  h  at  another,  24  and  8  splashes  were  heard  and 
spotlighted  within  5  m  of  the  boat.  Each  splash  was 
caused  by  adult  American  shad  in  a  tight  group  of 
2-6  individuals.  Observed  behavior  included  dash¬ 
ing  or  darting  (often  parallel  to  one  or  more  other 
fish),  tight  circling  with  another  fish  near  the  sur¬ 
face,  and  rolling  near  the  surface  parallel  to  another 
fish.  Gamete  release  was  not  observed. 

Spawning  activity  was  greatest  between 
30  April  and  17  June  and  between  2100  and  0100  h 
(Figs.  2  and  3).  In  PCA,  component  vectors  associ¬ 
ated  with  the  greatest  number  of  spawning 
splashes  (Fig.  4)  were  water  temperature  (positive 
relation  [+]),  DO  (negative  [-]),  and  current  velocity 
(-).  Regression  analysis  showed  significant  rela¬ 
tions  with  water  temperature  (+)  and  DO  (-),  cor¬ 
roborating  PCA  results  except  for  current  velocity 
(Table  4;  Figs.  5-8).  The  run  habitat  ranked  signifi¬ 
cantly  higher  in  spawning  activity  than  other  habi¬ 
tat  types;  pool  and  riffle  pool  were  lowest  (Table  5). 

Eggs  and  Larvae 

Postlarval  density  analyses  showed  no  evidence 
of  differential  catchability  by  gear  type.  Though 
the  mean  density  in  channels  was  2-3  times 
higher  with  bongo  than  plankton  seine  (Table  6), 


plankton-seine  data  came  from  a  single  location 
versus  three  sites  for  the  bongo.  In  addition,  habi¬ 
tat  density  ranks  were  the  same  regardless  of  gear 
used  to  sample  these  habitats  (Fig  5.).  Even  more 
conclusive  was  the  ratio  of  drift-net  to  bongo-net 
densities  for  prolarvae  and  postlarvae  (Table  7). 
Postlarvae  were  an  order  of  magnitude  denser  in 
drift  nets  than  bongo  nets  set  in  riffle  zones,  a 
result  opposite  of  expectations  were  net  avoidance 
by  postlarvae  possible  (drift  nets  have  much 
smaller  openings  than  bongo  nets).  The  greater 
numbers  of  postlarvae  likely  result  from  their 
occupying  a  higher  position  in  the  water  column 
(eggs  and  prolarvae  drift  nearer  to  the  substra¬ 
tum).  Possibility  of  greater  postlarval  escape  from 
the  benthic  sled  is  also  likely  owing  to  the  slow  tow 
velocities  needed  to  keep  the  sled  on  the  bottom. 
However,  sled  volumes  accounted  for  only  2%  of 
the  total  river  volume  sampled  (Table  3). 

American  shad  eggs  were  found  in  greatest 
densities  from  8  May  to  18  June  over  the  3- 
year  sampling  period,  prolarvae  8  to  17  June, 
and  postlarvae  14  May  to  9  July  (Fig.  9).  Over 
the  diel  cycle,  eggs  became  numerous  between 
2000  and  2400  h,  prolarvae  between  1800  and 
0030  h,  and  postlarvae  over  all  hours  sampled 
(0930-0230  h;  Fig.  10).  In  PCA,  American 
shad  eggs  showed  association  with  higher  cur¬ 
rent  velocities  and  shallower  river  depths 
than  larvae  (Fig.  11).  Riffle  habitat  was  more 
closely  and  slope  habitat  less  closely  aligned 


Fig.  2.  Relation  between  mean 
number  of  spawning  splashes 
per  5  min  and  date  for  adult 
American  shad  in  the  upper 
Delaware  River. 


Day  of  year  and  date 
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Fig.  3.  Relation  between  mean 
number  of  spawning  splashes 
per  5  min  and  time  of  day 
(Eastern  Daylight  Saving 
Time  [DST])  for  adult  Ameri¬ 
can  shad  in  the  upper  Dela¬ 
ware  River. 


Fig.  4.  Principal  component 
analysis  of  the  number  of 
adult  spawning  splashes 
(low,  moderate,  or  high)  in  re¬ 
lation  to  four  physical  habitat 
variables  in  the  upper  Dela¬ 
ware  River.  Mean  values  for 
low,  moderate,  and  high  num¬ 
bers  of  spawning  splashes  are 
plotted  as  L,  M,  and  H.  Vec¬ 
tors  are  drawn  to  twice  their 
calculated  length  (eigenvec¬ 
tor  output)  for  visual  en¬ 
hancement. 


Table  4.  Linear  and  quadratic  regression  analyses  among  five  physical  habitat  variables  and  number  of 
splashes  by  spawning  adult  American  shad  heard  on  the  upper  Delaware  River, 


Habitat  variable 

Regression 

type 

Degrees  of 
freedom 

F  value 

P 

Water  temperature 

linear 

1, 128 

10.8 

0.001** 

0.08 

Quadratic 

2, 127 

9.0 

0.011* 

0.12 

Dissolved  oxygen 

Linear 

1,  113 

7.6 

0.007** 

0.06 

Quadratic 

2, 112 

5.1 

0.12 

0.08 

Current  velocity 

Linear 

1, 128 

0.0 

0.90 

0.00 

Quadratic 

2, 127 

0.95 

0,17 

0.01 

Depth 

Linear 

1, 128 

0.12 

0.74 

0.00 

Quadratic 

2,127 

0.25 

0.54 

0.00 

Turbidity 

Linear 

1, 106 

3.70 

0.06 

0.03 

Quadratic 

2,105 

2.13 

0.12 

0.04 

•*  Significant  at  0.05  and  0.01  levels. 
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Fi§^.  5.  Habitat  suitability  index 
(HSI)  model  for  spawning 
adult  American  shad  (Stier 
and  Crance  1985)  superim¬ 
posed  on  relation  between 
mean  number  of  spawning 
splashes  per  5-min  interval 
and  surface  water  tempera¬ 
ture  for  adult  American  shad 
in  the  upi)er  Delaware  River. 
A  comparison  of  lower  pa¬ 
rameter  limits  for  minimum 
(0)  and  maximum  (1)  suitabil¬ 
ity  between  plotted  data  and 
the  Stier  and  Crance  (1985) 
model  is  presented  (inset). 


Fig.  6.  HSI  model  for  spawning 
adult  American  shad  (Stier 
and  Crance  1985)  superim¬ 
posed  on  the  relation  between 
mean  number  of  spawning 
splashes  per  5-min  interval 
and  current  velocity  for  adult 
American  shad  observed  in 
the  upper  Delaware  River.  A 
comparison  of  lower  parame¬ 
ter  limits  for  minimum  (0) 
and  maximum  (1)  suitability 
between  plotted  data  and  the 
Stier  and  Crance  (1985) 
model  is  presented  (inset). 
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Fig.  7.  Relation  between  mean 
number  of  spawning  splashes 
per  5-min  interval  and  river 
depth  for  adult  American 
shad  in  the  upper  Delaware 
River. 


Fig.  8.  Relation  between  mean 
number  of  spawning  splashes 
per  5-min  interval  and  tur¬ 
bidity  for  adult  American 
shad  in  the  upper  Delaware 
River. 
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Table  5.  Duncan*s  multiple  range  test  for  the  effect 
of  habitat  type  on  spawning  activity  by  American 
shad  in  the  upper  Delaware  River.  Error  degrees 
of  freedom  =  137,  F  =  5.75,  P<  0.001. 


Habitat  type 

n 

Mean  no.  spawning 
splashes 

Duncan 

group 

Run 

33 

23.1 

A 

Channel 

34 

10.5 

B 

SAV  shallows 

43 

10.2 

B 

Riffle  pool 

14 

5.1 

B 

Pool 

18 

1.6 

B 

with  egg  distribution  than  with  larval  distribu¬ 
tion.  Prolarvae  were  associated  with  river  depth 
(+)  and  current  velocity  (-)  to  a  greater  degree 
than  postlarvae  or  eggs.  Channel  and  pool  habi¬ 
tats  represented  prolarval  distribution  better 


than  postlarvae  or  eggs.  Postlarvae  were  associ¬ 
ated  with  higher  water  temperatures  and  lower 
turbidities  than  either  prolarvae  or  eggs.  Little 
difference  in  habitat  use  was  observed. 

Regression  analyses  corroborated  some  PCA- 
identified  relations  but  indicated  the  inverse  or  no 
relation  for  others  (Table  8;  Figs.  12-16).  Egg  den- 
sity  for  example,  showed  no  relation  to  current 
velocity  or  river  depth,  but  a  positive  relation  to 
sample  depth.  Prolarval  density  increased  only 
with  sample  depth  and  showed  no  relation  to  cur¬ 
rent  velocity.  Postlarvae  showed  positive  relations 
with  sample  depth  and  ambient  water  temperature 
among  the  six  variables  examined. 

Other  notable  findings  concern  the  range  over 
which  eggs  or  larvae  were  found  for  each  habitat 
variable  (Table  9).  Proleirvae,  for  example,  were 
found  over  a  much  shorter  portion  of  the  sampling 
season  than  eggs  or  postlarvae.  Eggs  were  found 
even  3  weeks  later  than  prolarvae.  Eggs  and  pro¬ 
larvae  escaped  capture  during  most  of  the  lighted 


Table  6.  ANOVA  and  Duncan  multiple  range  tests  for  effect  of  habitat  type  in  the  upper  Delaware  River 
on  American  shad  density  at  three  life  stages.  Postlarval  densities  are  also  partitioned  by  gear  type  to 
assess  potential  differences  in  gear  catchability  for  this  life  stage. 


Life  stage 
or  gear 

DF 

F 

P 

Duncan  multiple  range  test 
Mean  n  Habitat 

Group 

Model 

Error 

Egg 

4 

216 

0.8 

0.53 

0.396 

11 

Slope 

A 

0.213 

81 

Riffle 

A 

0.103 

103 

Channel 

A 

0.056 

16 

Pool 

A 

0.009 

10 

Riffle  pool 

A 

Prolarva 

4 

216 

1.1 

0.37 

0.0013 

103 

Channel 

A 

0.0008 

16 

Pool 

A 

0.0006 

10 

Riffle  pool 

A 

0.0004 

81 

Riffle 

A 

0.0000 

11 

Slope 

A 

Postlarva 

4 

216 

3.4 

0.01 

0.011 

10 

Riffle  pool 

A 

0.002 

81 

Riffle 

B 

0.002 

103 

Channel 

B 

0.001 

16 

Pool 

B 

0.000 

11 

Slope 

B 

Bongo  net 

2 

73 

0.5 

0.59 

0.0016 

62 

Channel 

A 

0.0007 

10 

Pool 

A 

0.0000 

4 

Slope 

A 

Plankton 

1 

26 

3.72 

0.06 

0.0106 

10 

Riffle  pool 

A 

seine 

0.0006 

18 

Channel 

A 

Drift  net 

0 

80 

0.0021 

81 

Riffle 

Benthic 

2 

33 

0.51 

0.61 

0.004 

23 

Channel 

A 

sled 

0.002 

6 

Pool 

A 

0.000 

7 

Slope 

A 
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Table  7,  Comparison  of  egg,  prolarval,  and  postlarval  densities  obtained  by  drift  or  bongo  net  in  riffle 

habitats  in  the  upper  Delaware  River.  _ 


Gear 

Date 

Eggs 

Prolarvae 

Postlarvae 

Drift  net 

15Jun92 

0.175 

0.005 

0.009 

0.082 

0.000 

0.000 

0.160 

0.004 

0.004 

16Jun92 

0.036 

0.000 

0.000 

0.010 

0.000 

0.010 

0.032 

0.000 

0.005 

17  June  92 

0.434 

0.000 

0.030 

0.876 

0.000 

0.038 

Mean 

0.226 

0.001 

0.012 

Bongo  net 

15  Jun  92 

0.349 

0.001 

0.000 

0.695 

0.000 

0.002 

16  Jun  92 

0.016 

0.000 

0.001 

0.054 

0.006 

0.006 

17  June  92 

1.306 

0.004 

0,006 

Mean 

0.484 

0.002 

0.003 

Ratio  drift  to  bongo 

0.47 

0.50 

4.0 

Fig.  9.  Relations  between 
American  shad  egg  and  larval 
density  and  date  on  the  upper 
Delaware  River. 
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Fig.  10.  Relations  between 
American  shad  egg  and  larval 
density  and  time  of  day  (DST^ 
on  the  upper  Delaware  River. 
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Fig.  11.  Principal  component 
analysis  for  the  densities  of 
American  shad  eggs,  prolar- 
vae,  and  postlarvae  as  they 
relate  to  six  physical  habitat 
variables  and  four  habitat 
types  in  the  upper  Delaware 
River.  Vectors  are  drawn  to 
twice  their  calculated  length 
(eigenvector  output)  for  vis¬ 
ual  enhancement. 
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Table  8.  Linear  regression  of  American  shad  egg 
and  larval  densities  on  six  physical  habitat 
variables  in  the  upper  Delaware  River.  DF  refers 


to  error  degrees  of  freedom. 


Habitat 

Life  stage 

variables 

Eggs 

Prolarvae 

Postlarvae 

Sample  depth 

DF 

227 

227 

227 

0.03 

0.03 

0.05 

Slope 

+ 

+ 

+ 

P 

0.007** 

0.009** 

0.0006*** 

River  depth 

228 

0.01 

228 

0.01 

228 

0.00 

Slope 

- 

+ 

+ 

p 

0.19 

0.26 

0.45 

Temperature 

DF 

227 

227 

227 

0.00 

0.01 

0.04 

Slope 

+ 

+ 

+ 

P 

0.36 

0.12 

0.003** 

Dissolved  oxygen 

DF 

226 

226 

226 

R^ 

0.00 

0.00 

0.01 

Slope 

+ 

- 

- 

P 

0.81 

0.21 

0.16 

Current  velocity 

DF 

222 

222 

222 

0.01 

0.01 

0.02 

Slope 

_ 

- 

- 

P 

0.13 

0.13 

0.06 

Turbidity 

DF 

191 

191 

191 

R^ 

0.00 

0.00 

0.00 

Slope 

- 

- 

“ 

P 

0.82 

0.77 

0.79 

**,  ** ‘Significant  at  0.01  and  0.001  levels. 


portion  of  the  diel  cycle.  Most  river  depths  were 
inhabited  by  both  eggs  and  larvae.  Larvae  were  not 
found  until  water  temperatures  reached  13°  C,  but 
eggs  and  larvae  were  found  up  to  the  highest  tem¬ 
peratures  sampled  (26.6°  C).  Eggs  and  larvae  were 
found  at  all  but  the  highest  (>1.1  m/s)  current 
velocities.  Egg  distribution  was  limited  to  turbidi¬ 
ties  <2.7  ntu,  while  larvae  were  restricted  to  tur¬ 
bidities  <2.1  ntu. 

Analysis  of  variance  for  the  effect  of  habitat 
type  (five  classes)  on  egg  and  larval  density 
showed  no  effect  on  eggs  or  prolarvae  but  a  signifi¬ 
cant  effect  on  postlarvae,  with  the  riffle  pool  habi¬ 
tat  showing  higher  densities  than  all  others  (Ta¬ 
ble  6).  The  Duncan  multiple  range  test  further 
showed  that  eggs  concentrated  in  slope  >  riffle  > 


channel  >  pool  >  riffle  pool;  prolarvae  in  channel 
>  pool  >  riffle  pool  >  riffle  >  slope;  and  postlarvae 
in  riffle  pool  »  riffle  >  channel  >  pool  >  slope. 

Juveniles 

Catchability  of  juvenile  American  shad  in 
SAV  habitats  ranged  from  52  to  80%;  the  aver¬ 
age  was  68%  (Table  10).  There  was  no  correla¬ 
tion  between  seine -haul  counts  and  dive  counts 
overall,  but  correlations  were  found  in  two  spe¬ 
cific  habitat  types:  SAV  (strongly  correlated)  and 
riffle  (Table  11). 

Juveniles  were  found  from  20  June  to  3  October 
and  at  all  hours  sampled  from  0353  to  2345  h 
(Figs.  17  and  18).  Principal  component  analysis  of 
the  1991  data  alone  showed  the  relative  abun¬ 
dance  of  juveniles  to  be  related  positively  to  water 
temperature  and  river  depth  but  negatively  to 
turbidity;  however,  these  relations  were  negated 
when  1992  data  (with  much  cooler  summer  water 
temperatures)  were  included  (Fig.  19).  Closer  re¬ 
lations  (positive  or  negative)  to  riffle  pool,  eddy, 
and  SAV  habitats  were  observed  than  to  riffle  and 
channel  habitats. 

Regression  analyses  with  dive  counts  alone  or 
dive  and  seine-haul  counts  showed  significant 
relations  for  two  of  six  physical  habitat  vari¬ 
ables  (Table  12;  Figs.  20-22).  Relative  abun¬ 
dance  of  juvenile  American  shad  was  directly 
related  to  water  temperature  and  inversely  to 
DO.  A  different  pattern  of  significant  relations 
emerged,  however,  when  regressions  were  per¬ 
formed  by  habitat  type  (Table  13).  Water  tem¬ 
perature  was  an  important  determinant  of  rela¬ 
tive  abundance  only  in  riffle  habitat,  while  river 
depth  became  important  in  SAV  habitat.  Cur¬ 
rent  velocity  and  turbidity  emerged  as  determi¬ 
nants  of  relative  abundance  almost  exclusively 
in  tributaries  (a  linear  relation  was  found  for 
current  velocity  in  riffle  pool).  Percent  SAV 
cover  was  important  in  eddy  (dive  counts  only) 
and  SAV  habitats. 

Analysis  of  variance  for  the  effect  of  habitat 
type  on  juvenile  relative  abundance  showed  no 
effect  (Table  14).  A  t-test  showed  no  difference 
between  the  first  and  second  (10  min  later)  dive 
count  on  grids  (regardless  of  habitat).  Means  ± 
standard  errors  were  33  ±  9  and  37  ±  11.  Neither 
means  (DF  =  48,  t  =  1.54,  P  =  0.13)  nor  variances 
(DF  =  48,  F  =  1.41,  P  =  0.24)  were  significantly 
different. 
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Yig,  12.  HSI  model  for  egg  and 
larval  American  shad  (Stier 
and  Crance  1986)  superim¬ 
posed  on  the  relations  be¬ 
tween  egg  and  larval  density 
and  surface  water  tempera¬ 
ture  observed  in  the  upper 
Delaware  River.  The  base  and 
peak  of  the  HSI  model  repre¬ 
sent  suitability  indices  of  0 
and  1,  respectively. 


Fig.  13.  Relations  between 
American  shad  egg  and  larval 
density  and  current  velocity 
in  the  upper  Delaware  River. 
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Fig.  16.  Relations  between 
American  shad  egg  and  larval 
density  and  tin-bidity  in  the 
upper  Delaware  River.  Fewer 
than  10  samples  were  taken 
at  higher  turbidity  readings 
(between  4.0  and  8.5  ntu); 
none  contained  eggs  or  lar¬ 
vae. 
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Table  9.  Comparison  of  the  measured  ranges  of  date,  time  of  day,  and  six  physical  habitat  variables  over 
which  American  shad  eggs  and  larvae  were  observed  and  sampled  in  the  upper  Delaware  River 


Variable  and 

Median 

Range 

life  stage 

observation 

Observed 

Sampled 

Date 

Eggs 

14  May 

16  Apr-09  Jul 

13  Apr-09  Jul 

97 

Prolarvae 

28  May 

02  May- 18  Jun 

13  Apr~09  Jul 

55 

Postlarvae 

16  Jun 

08  May-09  Jul 

13  Apr “09  Jul 

72 

Time  of  day  (h) 

Eggs 

2217 

1506-2628 

0954-2628 

67 

Prolarvae 

2214 

1543-2530 

0954-2628 

59 

Postlarvae 

2229 

0954-2628 

0954-2628 

100 

Sample  depth  (m) 

Eggs 

0.50 

0.20-4.5 

0-5 

86 

Prolarvae 

0.55 

0.25-4.0 

0-5 

75 

Postlarvae 

0.50 

0.25-4.8 

0-5 

91 

River  depth  (m) 

Eggs 

2.0 

0.2-4.8 

0-6 

77 

Prolarvae 

2.0 

0.2-4.8 

0-6 

77 

Postlarvae 

2.1 

0.5-5.0 

0-6 

80 

Temperature  (°  C) 

Eggs 

19.6 

8.2-26.6 

7.4-26.6 

96 

Prolarvae 

21.5 

13.0-26.2 

7.4-26.6 

69 

Postlarvae 

23.5 

13.0-26.6 

7.4-26.6 

69 

Dissolved  oxygen  (mg  /  L) 

Eggs 

9.5 

7.6-11.8 

7.6-14.1 

65 

Prolarvae 

9.3 

7.6-11.8 

7.6-14.1 

65 

Postlarvae 

9.1 

7.6-11.8 

7.6-14.1 

65 

Current  velocity  (m/  s) 

Eggs 

0.45 

0,01-1.10 

0.01-1.28 

76 

Prolarvae 

0.39 

0.01-0.98 

0.01-1.28 

76 

Postlarvae 

0.40 

0.01-1.02 

0.01-1.28 

80 

Turbidity  (ntu) 

Eggs 

1.2 

0.20-2.72 

0.20-8.30 

31 

Prolarvae 

1.1 

0.28-2.09 

0.20-8.30 

22 

Postlarvae 

1.1 

0.20-1.85 

0.20-8.30 

20 

®  Units  observed/  units  sampled  x  100. 
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Table  10.  Catchability  of  juvenile  American  shad  by  seining  in  SAV  habitats  (shallows  having  submerged 

aquatic  vegetation)  in  the  upper  Delaware  River. 

No.  juveniles  captiired _ 

Inner  net  haul®  no.  Outer  net  Catchability^ 

Date  Site  12  3  haul  Total  (%) 


llJul  91 

Zane  Grey 

14 

3 

0 

10 

27 

52 

20  Aug  91 

Milford 

8 

2 

10 

80 

09  Sep  91 

Dingman's  Ferry 

7 

1 

2 

2 

12 

58 

11  Sep  91 

Zane  Grey 

8 

1 

1 

0 

10 

80 

Total 

37 

5 

3 

14 

59 

68 

®  Inner  seine  net  hauls  were  made  within  an  encircling  (outer)  block  net. 

^Defined  as  the  percentage  of  fish  captured  by  the  first  inner  net  haul  of  the  total  fish  captured  by  all  net  hauls. 

Discussion 

Spawning  Adults 

Determination  of  habitat  requirements  for 
spawning  adult  American  shad  was  limited  in  this 
study  to  measurement  of  only  five  physical  habitat 
variables.  Habitat  use  could  have  been  quantified 
in  other  ways,  as  by  seine  deployment  and  elec¬ 
trofishing  techniques.  Such  methods  are  costly  in 
terms  of  time  and  resources,  however,  and  still 
would  not  ascertain  the  number  of  fish  actually 
spawning  in  a  specific  location.  Splashes,  on  the 
other  hand,  are  documented  indications  of  spawn¬ 
ing  activity  (Marcy  1972;  Chittenden  1976)  and  are 
readily  obtained  with  a  minimum  expenditure  of 
time  and  resources.  All  32  spotlighted  splashes 
were  verified  to  be  from  adult  American  shad  ex¬ 
hibiting  courtship  and  spawning  behavior.  The 
technique  was  thus  valid  and  useful  in  quantifying 
spawning  activity  in  this  species. 

A  major  problem  with  either  technique  is  scale 
and  habitat  resolution.  American  shad  seemed  to 


Table  11.  Pearson  correlation  coefficient  tests, 
overall  and  by  habitat  type,  for  paired  seine  and 
dive  counts  of  juvenile  American  shad  in  the 
upper  Delaware  River. 


Habitat  type 

n 

R 

P 

All 

39 

0.05 

0.78 

Submerged  aquatic 

18 

0.60 

0.009** 

vegetation 

Riffle 

5 

0.94 

0.02* 

Riffle  pool 

6 

-0.40 

0.43 

Eddy 

6 

0.32 

0.54 

Channel 

2 

Insufficient  data 

Tributary 

2 

Insufficient  data 

*,  **Significant  at  0.05  and  0.01  levels. 


spawn  over  large  areas,  both  longitudinally 
(>100  m)  and  laterally,  often  encompassing  several 
habitat  types.  Physical  habitat  variables,  espe¬ 
cially  depth  and  ctirrent  velocity,  vary  longitudi¬ 
nally  and  (especially)  laterally  from  a  single  loca¬ 
tion  on  shore.  However,  within  the  habitat 
limitations  acknowledged,  the  results  may  be  in¬ 
terpreted  reliably. 

Principal  component  analysis  identified  the  fol¬ 
lowing  parameters  as  important  correlates  of 
spawning:  water  temperature  (+)  >  DO  (-)  >  current 
velocity  (-).  Regression  analysis  identified  signifi¬ 
cant  relations  with  water  temperature  and  DO 
level.  Under  equilibrium  saturation  conditions  DO 
is  an  inverse  function  of  water  temperature.  Such 
oxygen-saturated  conditions  are  likely  to  exist  in 
the  river  during  the  spawning  season  (April-Jime) 
when  primaiy  productivity  (and  thus  photosynthe¬ 
sis)  has  not  yet  peaked.  Thus  a  water  temperature 
relation  can  be  inferred  from  the  DO  relation.  Dis¬ 
solved  oxygen  is  not  thought  to  affect  American 
shad  behavior  or  physiology  if  it  remains  above 
5  mg/L  (Chittenden  1969).  Accordiugly,  we  may 
ignore  DO  as  an  independent  variable  important  to 
this  life  stage. 

These  results  indicate  that,  other  than  tempo¬ 
ral  variables  (date  and  time  of  day),  water  tem¬ 
perature  is  the  primary  environmental  determi¬ 
nant  of  adult  spawning  activity  among  those 
variables  measured.  There  is  considerable  evi¬ 
dence  that  ambient  water  temperature  is  a  causa¬ 
tive  factor  in  the  initiation  of  migration  into  natal 
rivers  as  well  (Leggett  and  Whitney  1972). 
Though  not  demonstrated  specifically  in  Ameri¬ 
can  shad,  temperature  also  plays  a  role  in  stimu¬ 
lating  ovarian  maturation  in  fishes  (e.g.,  Sun- 
dararaj  and  Vasal  1976).  Our  observed  increase  in 
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Fig.  17.  Relation  between  juve¬ 
nile  American  shad  relative 
abundance  and  date  on  the 
upper  Delaware  River. 
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Fig.  18.  Relation  between  juve¬ 
nile  American  shad  relative 
abundance  8ind  time  of  day 
(DST^  on  the  upper  Delaware 
River. 
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Fig.  19.  Principal  component 
analysis  for  American  shad 
juvenile  abundance  in  rela¬ 
tion  to  five  physical  habitat 
variables  and  five  habitat 
tyi>es  in  the  upper  Delaware 
River.  Vectors  are  drawn  to 
twice  their  calculated  length 
(eigenvector  output)  for  vis¬ 
ual  enhancement. 


Table  12.  Linear  and  quadratic  regression  analyses  of  the  effect  of  six  physical  habitat  variables  on 
abundance  of  juvenile  American  shad  in  the  upper  Delaware  River  measured  by  both  dive  and  seine 
counts  or  by  dive  counts  only.  DF  refers  to  error  degrees  of  freedom. 

Both  dive  and  seine 

Dive  only 

Physical  variable 

Linear 

Quadratic 

Linear 

Quadratic 

River  depth 

245 

244 

204 

203 

0.09 

0.02 

0.01 

0.01 

Slope 

+ 

- 

-1- 

+ 

P 

0.15 

0.10 

0.20 

0.26 

Temperature 

DF 

248 

247 

206 

205 

0.03 

0.03 

0.04 

0.05 

Slope 

-1- 

+ 

-I- 

+ 

P 

0.01* 

0.03* 

0.002** 

0.004** 

Dissolved  oxygen 

DF 

248 

247 

206 

205 

0.01 

0.02 

0.02 

0.02 

Slope 

- 

- 

“ 

“ 

P 

0.12 

0.13 

0.04* 

0.09 

Current  velocity 

DF 

243 

242 

202 

201 

0.00 

0.00 

0.01 

0.01 

Slope 

- 

+ 

+ 

- 

P 

0.83 

0.90 

0.22 

0.40 

Turbidity 

DF 

231 

230 

191 

190 

P^ 

0.02 

0.02 

0.02 

0.02 

Slope 

- 

_ 

- 

- 

P 

0.03* 

0.07 

0.07 

0.11 

Percent  submerged  aquatic 
vegetation 

DF 

239 

238 

201 

200 

p2 

0.01 

0.01 

0.01 

0.01 

Slope 

+ 

- 

-1- 

+ 

P 

0.16 

0.27 

0.24 

0.46 

',  **  Significant  at  0.05  and  0.01  levels. 
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Fig.  20.  Relation  between  juvenile 
American  shad  relative  abundance 
and  water  temperature  in  the  upper 
Delaware  River. 
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Fig.  21.  Relation  between  juvenile 
American  shad  relative  abundance 
and  river  depth  in  the  upper  Delaware 
River. 
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Fig.  22.  Relation  between  juvenile  £ 
American  shad  relative  abundance  ^  200  1 
and  turbidity  in  the  upper  Delaware  | 
River.  g 
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spawning  activity  after  1900  h  until  0100  h 
(Fig.  3)  supports  findings  in  numerous  other  stud¬ 
ies  (e.g.,  Massman  1952;  Chittenden  1969).  Cur¬ 
rent  velocity  (discussed  below),  river  depth  (great¬ 
est  spawning  activity  at  <1  m;  Fig.  7),  and 
turbidity  (greatest  spawning  activity  at  <2  ntu; 
Fig.  8)  may  also  affect  spawning. 

The  HSI  models  for  spawning  adult  American 
shad  include  two  physical  habitat  variables,  mean 
surface  water  temperature  and  mean  water  veloc¬ 
ity  during  spawning  season  (Stier  and  Crance 


1985).  Superimposition  of  the  water  temperature 
model  onto  our  data  resulted  in  a  good  fit  at  the 
lower  end  of  the  distribution  (Fig.  5);  that  is, 
nearly  all  the  data  fell  within  the  SI  curve.  The 
data  did  not  fit  the  model  as  well  at  the  upper  end 
of  the  distribution,  however.  We  recommend  ex¬ 
tending  the  maximum-suitability  limits  to  24.5°  C 
accordingly. 

The  results  of  regression  analysis  (Table  4; 
Fig.  6)  and  PCA  (Fig.  4)  did  not  agree  for  current 
velocity.  In  PCA  an  inverse  relation  to  spawning 


Table  13.  Linear  (L)  and  quadratic  (Q)  relations  detected  from  analysis  of  regression  of  six  physical 
habitat  variables  on  abundance  of  juvenile  American  shad  by  habitat  type  in  the  upper  Delaware  River. 
Relations  listed  are  significant  (P  <  0.05)  for  both  dive  and  seine-haul  counts  unless  otherwise  noted. 


Habitat  type 

Habitat  variable 

ChanneP 

Eddy 

SAV^  Riffle 

Riffle  pool 

Tributary 

Temperature 
Dissolved  oxygen^ 
River  depth 

Current  velocity 
Turbidity 

Percent  SAV 

L,^Q« 

L,Q 

L,Q 

L,Q 

L" 

L,Q 

L,Q 

® There  were  no  significant  relations  determined  for  either  dissolved  oxygen  or  channel  categories. 
^  S AV  =  shallows  having  submerged  aquatic  vegetation. 

^  Dive  counts  only. 
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Table  14.  Comparison  of  dive-count  abundance  of 
juvenile  American  shad  in  six  habitat  types  in  the 
upper  Delaware  River. 


Habitat  type 

No.  of 
counts® 

Mean 

Duncan 

b 

group 

Eddy 

17 

42 

A 

Riffle  pool 

46 

35 

A 

Channel 

13 

31 

A 

Riffle 

46 

22 

A 

Submerged  aquatic 

74 

21 

A 

vegetation 

Tributary 

12 

16 

A 

^Each  count  =  mean  of  three  counts  over  a  S-m -square  grid  of 
substrate  by  a  diver  equipped  with  mask,  fins,  snorkel,  and 


^  ANOVA  results:  model  DF  5,  error  DF  202,  F  =  0.55,  P  =  0.74. 


activity  was  suggested,  while  regression  showed 
no  relation  and  a  flat  distribution  to  nearly 
0.75  nv^s  (Fig.  6).  Since  American  shad  seemed 
to  spawn  as  frequently  in  slow-moving  water 
(0-0.3  m/  s)  as  in  moderate  current  velocities 
(0.3-0.75  m/ s),  these  data  essentially  invali¬ 
date  any  lower  suitability  limit  and  indicate  an 
upper  maximum  suitability  limit  of  0.7  m/  s  ver¬ 
sus  the  0.9  m  /  s  limit  of  Stier  and  Crance  (1985). 
Problems  with  the  current-velocity  model  include 
scale  and  the  possibility  that  single-point  meas¬ 
urements  of  current  velocity,  even  in  habitats  with 
slow  or  still  waters,  may  not  represent  the  entire 
current-velocity  profile  experienced  by  fish 
spawning  there. 

Though  some  spawning  was  observed  in  all 
habitat  types  examined,  highest  activity  in  runs 
and  lowest  in  pools  and  riffle  pools  suggest  some 
habitat  selection  by  adult  spawners.  Physical  at¬ 
tributes  present  in  runs,  partially  characteristic 
of  channels  and  SAV  shallows,  and  absent  from 
riffle  pools  and  pools  may  explain  such  habitat 
preference.  A  combination  of  physical  charac¬ 
teristics  that  seems  to  be  avoided  by  spawning 
adults  is  slow  current  and  greater  depth.  One  or 
the  other  of  these  characteristics  is  present  in  the 
middle-ranked  habitat  types,  but  not  both  as  a 
general  observation.  However,  lowest-ranked 
habitat  types  are  characterized  by  both  deep  and 
slow  water.  Perhaps  either  swift  or  shallow  water 
(and  especially  both)  confers  higher  survival  to 
newly  spawned  eggs.  Alternatively  the  observed 
habitat  choice  may  be  unrelated  to  physical  vari¬ 
ables  and  reflect  some  unknown  ecological  selec¬ 


tive  pressure,  such  as  fewer  egg  predators  in  such 
habitats. 

Eggs  and  Larvae 

In  interpreting  the  results  of  egg  and  larval 
physical-habitat  data  analysis,  and  in  accordance 
with  Rexstad  et  al.  (1988),  we  place  confidence  in 
the  methods  of  analysis  as  follows:  regression 
analysis  >  PCA  >  range  analysis.  However,  PCA- 
identified  variables  were  rejected  if  regressions 
conflicted  with  the  type  of  relation  suggested  by 
PCA.  For  habitat-type  analysis,  we  place  full  con¬ 
fidence  in  the  ANOVA  and  the  further  discrimina¬ 
tion  provided  in  the  Duncan  multiple  range  test. 
With  these  criteria,  none  of  the  physical  habitat 
variables  (exclusive  of  temporal  factors  and  possi¬ 
bly  sample  depth)  seem  to  be  important  correlates 
or  determinants  of  American  shad  egg  density. 
There  was  no  significant  effect  of  habitat  type, 
which  coincides  with  the  nonmotile  nature  of  the 
passively  drifting  eggs. 

Prolarvae  related  strongly  only  to  sample  depth 
among  physical  variables  and  showed  no  ability  to 
select  habitat  type.  Postlarvae  related  strongly  to 
sample  and  river  depth  as  well  as  water  tempera¬ 
ture.  They  concentrated  in  riffle  pools  over  other 
habitat  types.  Riffle  pools,  characterized  by  mod¬ 
erate  depths  and  current  variable  in  velocity  and 
direction,  may  represent  a  refuge  from  drift  or  a 
better  opportunity  for  larvae  to  forage  and  feed. 
Since  sample  depth  was  significant  for  all  three 
life  stages  in  regression,  it  suggests  preferential 
habitat  use  of  the  deeper  river  zones  regardless  of 
river  depth.  Only  water  temperature  for  postlar¬ 
vae  stands  out  as  an  important  remaining  physi¬ 
cal  factor  identified  by  both  PCA  and  regression 
analysis. 

The  Stier  and  Crance  (1985)  HSI  model  for 
egg  and  larval  stages  of  American  shad  includes 
only  a  single  habitat  variable,  mean  surface 
water  temperature  during  the  period  when  eggs 
and  larvae  are  present.  A  minimum  DO  level  of 
5  mg/L  (easily  met  in  our  study  sites)  and  a 
maximum  mean  sediment  concentration  (mid¬ 
water)  of  100  ppm  (not  measured,  but  unlikely 
to  be  exceeded  at  any  study  location)  were  also 
specified.  For  evaluation  of  the  model,  egg  and 
larval  data  from  our  study  (only  positive  sam¬ 
ples,  i.e.,  where  density  >0)  were  superimposed 
on  the  SI  graph  for  the  same  stages  (Fig.  12). 
This  comparison  shows  that,  though  some  eggs 
were  found  below  the  thermal  limit  of  the  model 
(10°  C)  and  at  less  than  maximum  suitability  at 
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the  upper  thermal  limits  (26-27°  C),  moderate  to 
high  densities  of  eggs  fell  within  the  range  of  the 
model  for  maximum  suitability  (15-25°  C).  Pro- 
larvae  also  fell  within  the  range  of  maximum 
suitability  except  at  the  upper  end,  where  mod¬ 
erate  to  high  densities  were  found  at  >26°  C  (SI 
about  0,75).  Postlarvae  showed  no  reduction 
in  density  at  the  upper  end  of  the  thermal 
sampling  range  (26-27°  C).  These  data  indi¬ 
cate  that  the  model  needs  to  be  modified  for 
postlarvae,  and  possibly  for  prolarvae  as  well,  to 
extend  the  suitability  range  to  higher  tempera¬ 
tures.  Unfortunately,  our  data  do  not  allow 
specification  of  the  extent  of  change  needed. 

Our  data  also  point  to  the  importance  of  date  in 
determining  habitat  suitability  for  egg  and  larval 
stages.  This  finding  constitutes  a  seasonal  effect 
related  to  other  geographic  variables  such  as  lati¬ 
tude  and  river  under  study,  as  well  as  other  habi¬ 
tat  variables  such  as  ambient  water  temperature. 
Our  data  (Fig.  9)  indicate  maximum  suitability  of 
8  May  to  18  June  for  eggs,  8  May  to  17  June  for 
prolarvae,  and  14  May  to  9  July  for  postlarvae  on 
the  upper  Delaware  River.  Such  suitability  "win¬ 
dows"  are  clearly  specific  to  river  basins,  even  to 
river  segments.  The  window  limits  are  also  sub¬ 
ject  to  error  because  only  three  field  seasons  were 
used  to  determine  them. 

Time-of-day  distributions  for  eggs  and  larvae 
were  also  noteworthy  (Fig.  10).  The  distribution 
observed  for  eggs  may  result  from  a  combination 
of  three  factors:  (1)  time  of  spawning,  (2)  specific 
gravity  in  water  (slightly  heavier  than  water),  and 
(3)  tendency  of  currents  to  carry  suspended  eggs 
for  some  distance  (Massman  1952).  The  prolarval 
distribution  could  be  explained  by  a  vertical  swim¬ 
ming  response  to  light  (-),  although  sampling  error 
alone  might  be  responsible.  The  vertical  distribu¬ 
tions  depend  greatly  on  the  type  of  gear  used  to 
sample  each  life  stage.  Substrate  sampling  was  not 
accomplished  by  owe  gear,  for  example. 

The  distributions  of  egg  and  larval  stages  with 
regard  to  current  velocity,  sample  and  river  depth, 
and  tiirbidity  may  also  be  described  in  terms  of 
suitability  indices.  Such  interpretations  should 
not  be  construed  as  habitat  requirements,  how¬ 
ever,  because  no  causal  relations  have  been  dem¬ 
onstrated.  All  three  life  stages,  for  example, 
showed  "maximum  suitability"  ranges  in  current 
velocity  from  0.2  to  0.6  m/  s  (Fig.  13).  Our  data 
gave  little,  if  any,  indication  of  maximum  suitabil¬ 
ity  in  either  sample  or  river  depth,  however 
(Figs.  14  and  15).  Nearly  all  points  sampled  in  the 


water  column  at  all  river  depths  appeared  suitable 
for  all  three  life  stages.  Maximum  suitability  for 
turbidity  seemed  to  range  from  0.5  to  2.0  ntu 
(Fig.  16)  for  all  three  life  stages. 

Among  the  three  life  stages,  only  postlarvae 
showed  signific2Lntly  greater  habitat  use,  and  thus 
maximum  suitability,  for  a  particular  habitat 
type,  the  riffle  pool  habitat.  The  only  habitat  type 
sampled  where  no  eggs  or  larvae  were  obtained 
(zero  suitability)  was  the  shallow  inshore  SAV 
habitat. 

Juveniles 

Catchability  of  juvenile  American  shad  by 
seines  was  relatively  high  (52-80%)  in  SAV  habi¬ 
tats.  These  data  compare  to  a  mean  of  only  34%  in 
similar  habitat  observed  by  Parsley  et  al.  (1989). 
The  methods  used  to  determine  catchability  in 
SAV  could  not  be  used  in  other  habitats  because  of 
high  current  velocities,  deep  waters,  or  heteroge¬ 
neous  substrates.  Correlation  analysis  confirmed 
the  propriety  of  using  beach  seines  to  quantify 
juvenile  numbers  in  SAV  habitats  and  indicated 
its  justification  in  shallow  riffle  habitats  also 
(Table  11).  Because  of  relatively  nonturbid 
water  (90%  of  all  observations  at  <3  ntu)  with 
adequate  imderwater  visibility  (typically  3-5  m), 
underwater  visual  enumeration  of  juveniles  on 
established  quadrants  was  considered  more  reli¬ 
able  than  seine-haul  counts,  especially  in  the  case 
of  a  nonterritorial  schooling  species.  Follow-up 
dive  counts  (10  min  after  initial  counts)  were  con¬ 
sistent  with  initial  counts. 

The  results  of  regression  analyses  were  essen¬ 
tially  the  same  whether  dive  or  dive  and  seine  data 
were  used  (Table  12).  Since  PC  A  revealed  no  or 
only  weak  relations  to  either  physical  variables  or 
habitat  types,  conclusions  regarding  juvenile  habi¬ 
tat  use  are  based  primarily  on  the  regressions.  The 
following  habitat  variables,  in  order,  seemed  to  be 
important  overall  predictors  of  juvenile  abun¬ 
dance:  water  temperature  (+),  river  depth  (+),  tur¬ 
bidity  (-),  and  DO  (-). 

When  differentiated  by  habitat  type  (Table  13), 
however,  juvenile  abundance  was  affected  by  water 
temperature  only  in  riffles.  These  results  are  diffi¬ 
cult  to  explain  on  the  basis  of  energetics  because 
the  warmer  riffles,  where  juveniles  were  more 
abundant,  would  require  a  greater  expenditure  of 
energy  for  them  to  maintain  position  in  the  fast 
currents.  If,  however,  the  warmer  waters  hold  suf¬ 
ficiently  greater  concentrations  of  planktonic  or 
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drift  food  items,  the  energetic  benefits  may  out¬ 
weigh  the  costs. 

River  depth  affected  juvenile  abxmdance  pri¬ 
marily  in  SAV  habitat,  indicating  that  juveniles 
prefer  the  deeper  SAV  zones  somewhat  removed 
from  shore  and  shallows.  Turbidity  was  a  factor  (+) 
only  in  tributaries,  where  higher  turbidities  may 
have  been  associated  with  higher  concentration  of 
food  items.  Higher  current  velocities  associated 
with  the  same  habitat  type  may  also  have  been 
related  to  food  abundance.  Percent  SAV  cover  was 
found  important  only  in  SAV  and  eddy  habitats 
(though  SAV  was  present  to  lesser  degrees  in  other 
habitat  types),  indicating  that  aquatic  macro¬ 
phytes  provide  forage-related  advantage,  cover,  or 
some  other  resource  in  SAV  and  eddy  habitats  only. 
Since  no  effect  of  habitat  type  was  observed  on 
juvenile  relative  abundance  at  established  quad¬ 
rants,  juvenile  American  shad  seem  to  use  all 
habitat  tyi>es  studied  equally  to  their  advantage. 
Thus,  premigratory  juveniles,  in  contrast  to  earlier 
life  stages  and  spawning  adults,  may  be  habitat 
generalists. 

Proposed  HSI  Models  for 
Juveniles 

We  propose  three  HSI  models,  based  on  physi¬ 
cal  habitat  variables,  for  premigratory  juvenile 
American  shad  using  nursery  habitats  in  the 
Delaware  River.  The  models  as  proposed  are  of 
preliminary  nature  because  they  are  based  on 
only  two  seasons  of  observation.  They  are  also 
subject  to  modification  by  river  basin,  especially 
those  based  on  or  strongly  affected  by  temporal 
variables  (seasonal  and  time-zone  effects).  They 
are  based  on  underwater  dive  counts  only  because 
of  the  imcertainty  of  seine  efficiency  in  many 
habitat  types. 

Figures  20-22  show  juvenile  relative  abun¬ 
dance  at  all  study  sites  and  habitats  for  water 


temperature,  river  depth,  and  turbidity.  Dissolved 
oxygen  is  excluded  because  of  its  known  lower 
limit  for  juveniles  and  strong  relation  to  tempera¬ 
ture.  Maximum  (1)  and  minimum  (0)  SI  values, 
based  on  these  scatter  graphs,  are  listed  in  Ta¬ 
ble  15.  Because  most  of  the  input  variables  to 
these  models  were  shown  to  be  habitat  specific 
(Table  13),  consideration  should  be  given  to  the 
range  and  diversity  of  habitat  types  in  the  river 
system  imder  study  before  applying  the  models. 
Further  study  is  needed  to  determine  upper  and 
lower  limits  to  the  models  where  insufficient  data 
precluded  their  specification. 

The  HSI  models  for  American  shad  originally 
developed  by  Stier  and  Crance  (1985)  and  modified 
or  extended  to  include  additional  life  stages  by  this 
study  may  be  applied  broadly  to  all  river  systems 
with  caution.  Our  data  indicate,  however,  that 
within  a  specific  river  basin,  suitable  ranges  of 
physiceJ  factors  for  given  life  stages  may  be  either 
broader  or  narrower  than  previously  presumed.  We 
also  found  that  such  suitability  ranges  can  be  habi¬ 
tat  specific  within  river  segments.  These  changes 
may  be  due  to  differences  in  food  abimdance,  preda¬ 
tors,  or  other  unknown  or  unexamined  biotic  factors 
within  specific  habitats.  If  the  range  and  proportion 
of  habitat  types  in  each  river  system  of  interest  are 
known,  HSI  models  for  American  shad  may  be 
applied  cautiously  on  a  comparative  basis  among 
river  systems. 
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Table  15.  Maximum  (1)  and  minimum  (0)  suitability  indices  (SI)  for  juvenile  American  shad  in  nursery 
habitats  for  three  habitat  variables  identified  as  predictors  of  relative  abundance  in  the  upper  Delaware 
Riven  Missing  values  represent  undetermined  limits  because  of  insufficient  observations. 


Lower-limit  SI  Upper-limit  SI 


Habitat  variable 

0 

1 

1 

0 

Water  temperature  (°  C) 

17 

19.5 

24.5 

River  depth  (m) 

0.2 

1.2 

Turbidity  (ntu) 

0.0 

0.0 

2.2 

5.0 
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our  fish,  wildlife,  and  biological  diversity:  preserving  the  environmental  and  cultural 
values  of  our  national  parks  and  historical  places;  and  providing  for  the  enjoyment 
of  life  through  outdoor  recreation.  The  Department  assesses  our  energy  and  mineral 
resources  and  works  to  ensure  that  their  development  is  in  the  best  interests  of  all 
our  people  by  encouraging  stewardship  and  citizen  participation  in  their  care.  The 
Department  also  has  a  major  responsibility  for  American  Indian  reservation 
communities  and  for  people  who  live  in  island  territories  under  U.S.  administration. 


